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ABSTRACT 


This  report  presents  and  interprets  the  predictions  of  an  unsteady  aerodynamic  prediction 
method  known  as  the  Sonic  Box  method.  Illustrations  are  given  on  how  the  program  interprets 
input  modal  data,  the  program’s  techniques  for  smoothing  certain  input  and  output  data, 
convergence  of  the  numerical  results,  and  comparisons  of  predicted  results  with  experiments. 
It  is  shown  bow  the  present  program  requires  the  user  to  devote  some  care  to  defining  input 
mode  shapes,  however  this  problem  can  be  removed  by  a  simple  modification.  Generally 
speaking,  the  program’s  current  limit  of  fifty  boxes  in  any  one  direction  is  sufficient  to  obtain 
satisfactory  convergence,  with  the  exception  of  pitching  moments  on  cropped  deltas.  In  this 
respect  other  modifications  are  apparent  which  could  improve  convergence  at  the  cost  of 
increased  complexity.  Agreement  with  experiment  was  generally  qualitatively  good,  but 
illustrated  the  need  for  further  optimization  of  the  method  as  well  as  the  lack  of  experimental 
data  of  the  type  and  quality  desired  for  correlation. 
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SECTION  I 
INTRODUCTION 


The  development  of  a  numerical  method  for  calculating  the  aerodynamic  pressures  and 
generalized  forces  on  vibrating,  planar  lifting  surfaces  at  Mq,  =  1.0  is  found  in  Reference  1. 

The  solution  was  based  on  the  linearized  transonic  equations  of  Reference  2  which  were 
satisfied  by  a  distribution  oi  velocity  potential  doublets  over  the  wing  surface  and  its  wake. 
Generally  speaking  the  linearized  equation  holds  only  for  “high”  frequency  motions  of  a 
wing  with  *  ‘highly”  swept  leading  edges.  The  problem  was  further  simplified  by  restricting 
attention  to  wings  with  straight,  unswept  trailing  edges  which  then  allowed  the  wake  effects 
to  be  neglected. 

By  covering  the  wing  with  a  fine  grid  of  square  boxes  and  assuming  that  the  velocity  poten¬ 
tial,  <f> ,  was  constant  in  each  box,  the  integral  equation  for  the  velocity  potential  was  replaced 
by  a  system  of  linear  algebraic  equations  for  the  constant  value  in  each  box. 


f rhmrnrjnthnnirrrj 

SKETCH  OF  WINS  AH D  APPROXIMATE  6RIDV0RK 


Although  Reference  1  contained  a  complete  summary  of  the  development  of  the  equations 
and  a  resulting  FORI  RAN  IV  computer  program,  it  did  not  fully  demonstrate  the  application 
and  usefulness  of  the  computer  program.  The  purpose  of  this  report  is  to  accomplish  that 
demonstration. 
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SECTION  n 

INTERPRETATION  OF  MODE  SHAPE  INFORMATION 


The  primary  purpose  of  the  “Sonic  Box’  *  computer  program  is  to  generate  generalized 
forces  for  use  in  flutter  or  gust  response  calculations  for  surfaces  in  transonic  flow.  It  is 
very  important  that  the  user  understand  the  various  non-dimensionalizing  and  normalizing 
factors  that  are  in  the  program.  The  user  then  can  be  assured  that  the  mode  shapes  applied  to 
generate  the  aerodynamic  forces  are  identical  in  every  respect  to  those  used  to  generate 
generalized  mass  terms  for  flutter  calculations.  For  instance,  the  sonic  box  program  accepts 
mode  shape  information  in  either  of  two  ways: 

a.  Coefficients  of  modal  polynomials 

b.  Discrete  deflections  at  various  coordinate  points 

Questions  naturally  arise  in  the  user’s  mind  as: 

“Are  the  modal  polynomials  in  terms  of  dimensional,  non-dimensional  or  mixed  variables? 

When  I  put  in  a  deflection  of  1.0,  does  it  mean  one  unit  of  length  or  one  root  chord?” 

Appendix  I  presents  a  method  of  answering  these  questions  for  the  user  of  this  or  any  similar 
aerodynamic  computer  program  who  may  not  be  totally  familiar  with  the  program.  The  method 
of  Appendix  I  will  be  used  in  the  remainder  of  this  section  to  determine  exactly  the  forms 
required  for  the  proper  interpretation  of  modal  data. 

To  apply  the  method  a  dimensional  (ft,  meters,  and  so  forth)  coordinate  system  X,  y,  z  is 
adopted. 


A  computer  program  will  be  employed  which  uses  coefficients  &mn  of  a  polynomial  surface 

to  define  mode  shapes,  but  it  is  uncertain  whether  the  independent  or  dependent  variables  are 
divided  by  the  root  chord.  Therefore  the  equation  is  used  for  the  polynom  ial  surfuC8  in  tsrms 
of  two  unknowns  r  and  t. 


■—  -  I  Amn  <7  /bt)m  (7/b’)n  (i ) 

b  m,n 

For  instance,  if  the  powers  r  and  t  are  both  zero,  then  the  coefficients  should  be  deter¬ 
mined  by  using  the  dimensional  parameters  x,  y,  z.  On  the  other  hand,  if  r  and  t  are  both  1.0, 

then  A  should  be  determined  using  the  non-dimensional  forms  X/b,  y/b,  z/b.  It  should  «Js© 
inn 

be  pointed  out  that  it  may  be  possible  for  r  and  t  to  be  unequal  or  even  to  be  fractions. 
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The  generalized  forces  from  two  successive  computer  runs  on  the  same  wing  at  zero 

frequency  can  be  used  to  determine  r  and  t.  In  the  first  run  a  plunge  mode  is  defined  by  the 

coefficient  A  being  nonzero  with  all  other  A  =  0.0.  In  the  same  run  a  rotation  mode  is 
OO  °  mn 

also  defined  by  making  the  coefficient  A^q  nonzero  with  all  the  remaining  A^n  =  0.0.  In  a 
second  run  all  input,  is  kept  the  same,  except  the  dimensions  of  the  wing  are  scaled  by  a  factor 
of  Rg  =  b^Vb*2\  If  the  relationship  of  the  lifts  and  moments  is  made  known  between  the  two 
runs  by  Rl  =  L21(1)/L21(2)  andRM  =  L22(1)/L22(2).  it  then  turns  out  that  the  unknown  powers 

r  and  t  must  be: 


r 


1/2  (  1+2 


log  Rl 
log  Rg 


|pq  rm  \ 

lo  g  R  b  ' 


(2) 


log  Rl 
log  Rg 


lo0  rM 
log  Rg 


(3) 


In  deriving  the  above  forms  only  steady  state  motion  is  considered,  k  =  0.0.  It  turns  out, 
however,  that  the  Sonic  Box  computer  program  cannot  accept  a  frequency  of  zero  so  that  the 
following  runs  for  a  63.5  delta  wing  were  made  for  nonzero  (but  very  small)  frequencies.  Run 
number  1  was  made  with  w  =  0.0001  CPS.  Run  number  2,  with  a  wing  one  tenth  as  large,  was 
run  with  w  =  0.001  CPS  to  keep  k  the  same  for  both  runs.  Table  I  contains  the  important  input 
and  output  information.  Appendix  n  and  Reference  1  contain  explanations  of  the  input  format. 

Table  I  reveals  that  Rg  =  10.0,  RL  =  10.0,  and.  RM  =  100.0.  Therefore,  Equations  2  and 

3  give  r  =  0.5  and  t  =  0.  Then  when  the  user  furnishes  the  computer  program  with  coeffi¬ 
cients  Ajnn  the  generalized  forces  are  printed  out  as  if  the  user  is  actually  supplying  a 

polynomial  surface  of  the  form: 


z  r*  .  *vfn  -N-n 

7 =="  =  2.  Amn  x  y 

v  D  m,n 


where  b,  x,  y,  z  are  dimensional  quantities 


Having  answered  the  question  about  the  meaning  of  the  coefficients  A  when  they  constitute 

the  modal  input  data,  it  should  be  determined  how  the  computer  program  interprets  modal 
data  when  deflections  are  furnished  at  discrete  coordinate  points.  That  is,  certain  questions 
must  be  answered: 


"When  coordinates  are  put  in  should  they  be  dimensional  coordinates  X,  y,  or  the  non- 
dimensional  coordinates  x/b,  y/b?  When  a  number  is  put  in  for  the  deflection,  is  it  entered 
in  z  or  z/b?” 


To  answer  the  questions  the  technique  shown  in  Appendix  I  must  be  used  and  two  more 

computer  runs  made  which  are  similar  to  the  runs  made  to  determine  the  meaning  of  the 

coefficients  A _ . 

mn 


For  the  first  run  at  a  frequency  of  0.0001  CPS  b  =■  100.0  was  put  in  to  represent  the  root 
chord  and  1.0  to  represent  the  deflection  (plunge)  at  the  apex  (x  »  0.0)  and  at  some  other  x 
coordinate  (say  x  =  p).  In  the  same  run  a  rotation  mode  was  also  put  in  by  placing  a  deflection 
of  0.0  at  the  apex  and  1.0  at  .  For  a  second  run  at  a  frequency  of  0.001  CPS  and  b  -  10.0  a 
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deflection  of  1.0  was  first  put  in  at  the  apex  and  at  £  to  simulate  a  plunge  mode.  In  the  same 
run  a  deflection  of  0.0  was  then  put  in  at  the  apex  arid  1.0  at  @  to  simulate  a  rotation.  Table  II 
contains  the  important  input  and  output  information. 

Table  n  reveals  that  Rfi  -  10,  *  10,  and  R^  =  100.  Therefore  Equations  2  and  3  give 

r  =  0.5  and  t  =  0.0.  Then  when  the  user  furnishes  numbers  for  a  set  of  coordinates  and  de¬ 
flections,  the  generalized  forces  are  printed  out  as  if  he  is  actually  supplying  data  of  the  form: 

Coordinates:  x,  y  _ 

Deflections:  zf/vb 

where  b,  x,  y,  z  are  dimensional  quantities. 

It  is  important  to  note  that  similar  analyses  can  be  performed  to  show  that  pressures  are 

z 

printed  out  by  the  program  as  if  the  modal  data  were  furnished  in  the  form  ras  a  function 
of  x  and  y. 
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SECTION  m 

SMOOTHING  OF  INPUT  AND  OUTPUT 


The  Sonic  Box  computer  program  uses  smoothing  techniques  in  two  ways.  First,  when  modal 
deflections  are  given  by  the  user  in  terms  of  discrete  deflections  at  coordinate  points  the 
program  automatically  fits  a  least  squared  error  polynomial  surface  to  the  deflection  data. 
Secondly,  when  the  program  calculates  velocity  potentials  at  the  center  of  each  box  it  again 
automatically  fits  a  least  squared  error  polynomial  surface  to  the  data.  The  equation  of  the 
polynomial  surface  for  <f>  is  then  used  to  evaluate  the  unsteady  pressures  on  the  wing. 

Figures  1-4  illustrate  the  smoothing  of  deflection  data.  Figure  1  is  a  sketch  of  a  70°  delta 
wing  which  will  be  referred  to  extensively  in  later  sections  of  this  report  and  which  was 
subjected  to  various  rigid  body  and  flexible  mode  shapes.  Figure  2  illustrates  the  results  of 
applying  a  mode  shape  of  a  half  sine  wave  in  the  chordwise  direction.  Discrete  deflections 
were  specified  at  0,25,50, 75,  and  100  percent  of  the  root  chord.  The  least  squares  polynomial 
was  plotted  at  various  locations  and  can  not  be  distinguished  from  the  desired  half  sine  wave 
for  all  practical  purposes.  Figure  3  is  a  similar  plot  for  a  full  sine  wave  in  the  chordwise 
direction.  Here  the  least  square  fit  is  less  adequate  and  is  seen  to  distort  the  mode  shape 
slightly.  Figure  4  is  a  similar  plot  for  3/2  sine  waves  in  the  chordwise  direction  and  reveals 
a  progressively  less  accurate  approximation  with  proportionately  greater  number  of  deflection 
points  given.  Table  m  summarizes  the  important  input  data. 

Figures  5-15  illustrate  the  smoothing  of  calculated  velocity  potentials.  Figure  5  is  a  sketch 
of  an  aspect  ratio  1.5  delta  that  was  the  subject  of  calculation  for  modes  of  plunge  and 
rotation.  Figure  6  is  a  plot  of  the  chordwise  distribution  of  Im  (<f>)  for  a  spanwise  station  of 
1.0  ft  or  26  2/3  percent  of  the  semispan  for  a  plunging  motion.  Figure  6  is  drawn  for  suc¬ 
cessively  increasing  the  number  of  boxes  along  the  root  chord  from  5  to  45.  The  dots  on  the 
figure  represent  the  initially  calculated  values  of  the  velocity  potential,  whereas  the  lines 
represent  the  least  squares  polynomial  fits.  For  as  few  as  fifteen  boxes  along  the  root  chord, 
the  curve  fitting  procedure  gives  reasonable  predictions  for  d>.  Figure  7  is  a  similar  plot 
for  Real  (<j>)  in  rotational  motion.  The  same  progression  is  evident,  and  fifteen  boxes  along 
the  root  chord  appear  to  be  a  sufficient  number  to  get  reasonable  behavior  for  <£.  Figure  8 
is  merely  a  blown  up  presentation  of  the  bottom  curve  in  Figure  6  and  is  a  good  illustration 
of  the  behavior  of  near  the  leading  edge.  Note  that,  for  the  chordwise  line  selocted,  the 
leading  edge  is  at  26  2/3  percent  of  the  root  chord.  Figure  9  gives  an  isometric  view  of  the 
initially  calculated  values  for  Im  (<fi)  in  the  plunge  mode,  and  Figure  10  gives  an  isometric 
view  of  how  <f>  is  smoothed  by  the  program.  See  Table  IV  for  the  input  data. 

Figure  11  is  a  sketch  of  an  aspect  ratio  3  rectangular  wing  which  was  the  subject  of  analysis 
and  experiment  reported  in  Reference  3.  Figure  12  shows  the  first  bending  mode  which  has 
been  normalized  to  give  a  tip  deflection  of  4.098  for  reasons  of  convenience  which  will  be 
discussed  in  Section  V.  See  Table  V  for  the  input  data.  Figure  13  is  a  plot  of  die  chordwise 
variation  of  Im  (tp)  at  a  50  percent  semispan  station  for  increasing  numbers  of  boxes  along 
the  root  chord.  Figures  14  and  15  are  similar  plots  for  larger  numbers  of  boxes.  Note  that, 
for  eleven  or  more  boxes  along  the  root  chord,  the  velocity  potential  shows  erratic  behavior 
near  the  leading  edge  which  actually  grows  more  severe  as  the  number  of  boxes  is  increased. 
This  is  probably  due  to  the  basic  assumption  of  this  Sonic  Box  method  that  <f>  is  constant  over 

any  one  box,  whereas,  near  the  leading  edge  actually  should  approach  infinity.  A  correc¬ 
tion  factor  is  applied  for  rectangular  wings  in  the  existing  program  by  multiplying  <f>  in  the 
first  row  of  boxes  by  -J-  (Reference  1).  However,  apparently  this  correction  is  inadequate. 

A  revised  method  should  be  generated  which  allows  <f>  to  vary  within  any  box,  particularly 
near  the  leading  edge. 
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Figure  16  illustrates  the  smoothing  of  <f>  for  a  plunge  motion  of  the  70-degree  delta  wing  of 
Figure  1;  see  Table  Ql  for  a  summary  of  the  input  data. 
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SECTION  IV 
CONVERGENCE 


One  of  the  crucial  properties  of  any  numerical  method  is  that  the  answers  converge. 
Figures  17-53  reflect  an  exhaustive  study  of  the  convergence  of  the  Sonic  Box  program  for 
the  aspect  ratio  1.5  delta  wing  of  Figure  5,  the  aspect  ratio  3  rectangular  wing  of  Figure  11, 
the  aspect  ratio  2  delta  wing  of  Figure  27,  the  cropped  deltas  of  Figures  30  and  31,  md 
particularly  the  70 -degree  delta  of  Figure  1. 

In  Section  II  it  was  shown  that  when  the  user  employs  polynomial  coefficients  for  mode 

shapes  or  gives  coordinates  and  discrete  deflections,  the  generalized  forces  are  printed  out 
as  if  z4/d  is  a  function  of  x  and  y.  By  a  similar  analysis  (Appendix  I)  it  can  be  shown  that 
the  programs  print  out  pressures  as  if  the  user  wrote  z/b  as  a  function  of  X  and  % 

In  this  report  the  actual  pressure  coefficient  will  be  denoted  by  Ap/q  and  the  value  of 
“pressure”  printed  out  the  program  by  (Ap/q)*.  Similar  differentiations  with  respect  to  the 
generalized  force  coefficients  will  be  made. 

1.  THE  ASPECT  RATIO  1.5  DELTA 

The  wing  considered  is  shown  in  Figure  5.  The  important  input  data  is  summarized  in 
Table  IV.  Figure  17  is  a  plot  of  the  behavior  of  the  absolute  value  of  the  generalized  force 
coefficients  Lf.  as  the  number  of  boxes  along  the  root  chord  varied  from  5.  15,  25,  35.  and 

45.  In  this  report  the  notation  of  Reference  1  was  adopted: 


L'i  =  q 


where: 


ffs  dS 


Ap^  =  the  pressure  difference  -  pu  in  the  i  mode 

th 

z’j  =  the  deflection  in  the  j  mode 


Figure  17  shows  fairly  rapid  convergence  when  a*  least  fifteen  boxes  are  used  along  the  root 
chord  and  that  twenty  boxes  are  probably  sufficient  for  reasonable  accuracy.  Figure  18  is  a 
chordwise  plot  of  the  imaginary  part  of  the  pressure,  Im  (Ap/q)*,  at  a  spanwise  station  of 
26  2/3  percent  for  varying  numbers  of  boxes  along  the  root  chord.  Even  though  the  generalized 
forces  seemed  to  be  reasonably  well  converged  for  twenty  or  more  boxes,  there  is  considerable 
variation  of  the  pressure  near  the  trailing  edge  as  the  number  of  boxes  is  changed.  Figure  19 
is  a  similar  plot  for  Real  (Ap/q)*  for  the  rotational  mode  with  similar  results  for  the  trailing 
edge  pressures.  The  oscillatory  behavior  of  the  trailing  edge  pressures  leads  one  to  doubt 
the  true  convergence  so  that  a  second  computer  run  was  made,  varying  the  number  of  boxes 
in  increments  of  1.0.  Figure  20  illustrates  the  results  which  indicate  that  perhaps  as  many 
as  thirty  boxes  are  required  along  the  root  chord  for  convergence  of  L*^  and  that  the  oscilla¬ 
tory  behavior  of  the  trailing  edge  pressures  does  not  influence  the  generalized  forces  to  any 
great  extent.  The  second  computer  run  also  revealed  that  the  trailing  edge  pressures  can 
change  drastically  when  the  number  of  boxes  along  the  root  chord  is  changed  by  only  1.0. 
Figure  21  plots  1m  (Ap/q)*  in  the  plunge  mode  for  14,  15,  and  16  boxes  along  the  root  chord. 
Note  the  rapid  oscillation  in  (Ap/q)*.  The  results  in  Figure  22  are  similar  except  the  data 
are  identical  for  24  and  25  boxes  with  a  large  jump  in  the  predicted  trailing  edge  pressure 
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as  the  number  of  boxes  is  increased  to  26.  Figures  23  and  24  show  that  predicted  pressure 
finally  begins  to  settle  down  as  the  number  of  boxes  is  increased  to  36  and  then  46. 

2.  THE  ASPECT  RATIO  3  RECTANGLE 

Figure  25  is  a  plot  of  the  chordwise  pressure  distribution  over  the  aspect  ratio  3  rectangular 
wing  of  Figure  11  for  the  wing  bending  mode  of  Figure  12.  See  Table  V  for  input  data.  The 
Im  (Ap/q)*  is  plotted  at  a  50  percent  spanwise  station  for  6,  11,  15,  and  19  boxes  along  the 
chord.  The  pressures  in  Figure  25,  as  well  as  those  in  Figure  26  for  larger  numbers  of  boxes, 
appear  to  be  well  converged.  However,  for  later  reference,  note  the  apparently  spurious 
behavior  of  the  pressure  near  the  trailing  edge.  In  the  absence  of  numerical  problems,  one 
would  not  expect  this  method  to  predict  increasing  pressures  over  the  aft  10  percent  of  the 
chord.  Also  note  that  the  predicted  pressures  in  the  vicinity  of  10  -  30  percent  chord  are 
probably  too  low.  Other  theoretical  and  experimental  results  indicate  higher  pressures  near 
the  1/4  chord. 

3.  THE  ASPECT  RATIO  2  DELTA 

Figure  27  is  a  sketch  of  an  aspect  ratio  2.0  delta  wing  which  was  the  subject  of  unsteady 
1  essure  measurements  reported  in  Reference  4.  The  model  was  oscillated  in  rotation  about 
an  axis  at  approximately  71  percent  of  the  root  chord.  The  important  input  parameters  are 
summarized  in  Table  VI. 

Figure  28  is  a  plot  of  the  moment  due  to  rotation  versus  the  number  of  boxes  along  the 
root  chord.  Note  that  the  predicted  values  are  converging  to  a  final  answer,  but  not  partic¬ 
ularly  rapidly.  Figure  29  is  a  plot  of  the  chordwise  pressure  along  a  36  percent  semispan 
station  in  the  form  used  in  Reference  4.  There  is  a  slight  oscillation  in  the  pressures  with 
increasing  numbers  of  boxes.  Note  the  fairly  good  agreement  with  these  results  and  those 
from  the  Jones’  low  aspect  ratio  theory  taken  from  Reference  4. 

4.  THE  CROPPED  DELTAS 

Figures  30  and  31  are  sketches  of  two  ’‘cropped”  delta  wings  which  were  the  subject  of 
oscillatory  rotational  measurements  reported  in  Reference  5.  Each  of  the  wings  was  oscillated 
about  a  forward'  and  an  aft  axis,  as  shown  in  Figures  30  and  31.  The  significant  input  data  is 
summarized  in  Table  VII. 

Figure  32  presents  a  plot  of  the  coefficient  LJ^  (in  this  case,  moment  due  to  rotation) 

versus  the  number  of  boxes  along  the  root  chord  of  the  aspect  ratio  1.5  cropped  delta.  Al¬ 
though  the  data  for  the  aft  rotation  axis  is  somewhat  more  converged  than  for  the  forward 
rotation  axis,  it  is  apparent  that  forty  boxes  along  the  root  chord  is  not  sufficient  for  con¬ 
vergence  here.  The  same  comments  apply  for  the  aspect  ratio  2.0  wing  in  Figure  33,  the 
least  converged  being  the  real  part  of  the  moment  for  the  forward  rotation,  axis. 

C  rp TJ n  r\PT  m  A 

O.  *****  IV  UCiUl/1 


Extensive  convergence  checks  were  performed  for  the  70°  delta  wing  of  Figure  1,  using 
four  modes:  plunge,  rotation  about  the  apex,  a  half  sine  wave  flexible  mode  in  the  chordwise 
direction,  and  a  full  sine  wave  flexible  mode  in  the  chordwise  direction.  Table  HI  summarizes 
the  input  data. 

Figures  34-38  are  for  the  plunge  mode.  Figure  34  shows  the  behavior  of  the  lift  due  to 
plunge  versus  the  number  of  boxes  along  the  root  chord.  It  appears  that  thirty  or  more  boxes 
are  required  to  obtain  convergence  within  about  10  percent  of  the  final  answer.  To  show 
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that  the  velocity  potential  has  settled  down,  Figure  35  was  prepared  to  show  the  chordwise 
variation  of  near  the  root  chord  for  48,  49,  and  50  boxes  along  the  root  chord.  The  results 
for  the  three  cases  could  not  be  distinguished.  The  imaginary  part  of  the  pressure  near  the 
root  chord  is  shown  in  Figure  36.  The  results  are  certainly  converged  for  48,  49,  and  50  boxes, 
except  for  some  -light  oscillation  at  the  trailing  edge.  When  one  moves  outboard  to  50  percent 
span  the  pressure  results  in  Figure  37  are  still  well  converged  with  twenty  boxes  giving  at 
least  qualitative  agreement.  Figure  38  is  a  spanwise  plot  of  the  pressure  in  plunge,  taken  at 
x  =  2/3  of  the  root  chord.  Again,  the  convergence  for  48,  49,  and  50  boxes  and  the  qualitative 
agreement  for  twenty  boxes  are  apparent. 

Figures  39-44  are  for  the  rotation  mode.  Figure  39  is  a  plot  of  the  moment  due  to  rotation 
about  the  apex.  As  has  been  characteristic  so  far  of  rotation  modes,  the  moment  is  not 
particularly  well  converged.  The  program’s  upper  limit  of  fifty  boxes  does  not  give  assurance 
of  a  converged  solution.  Figure  40  illustrates  that,  even  though  the  integrations  which  result 
in  generalized  forces  have  not  converged  for  fifty  boxes,  the  predicted  values  of  the  velocity 
potential  near  the  root  are  reasonably  well  predicted  for  as  few  as  twenty  boxes.  In  fact,  the 
<f>  results  for  48,  49,  and  50  boxes  could  not  be  distinguished.  The  same  comments  apply  at 
the  50  percent  span  station  as  shown  in  Figure  41.  One  would  certainly  say  that  the  velocity 
potential  results  are  convergent.  The  pressure  results  near  the  root  are  shown  in  Figure  42 
to  be  fairly  well  converged  with  the  usual  oscillation  near  the  trailing  edge.  The  pressures 
at  50  percent  of  the  semispan  in  Figure  43  are  slightly  less  converged  but  still  acceptable. 
Again  the  results  for  as  few  as  twenty  boxes  along  the  root  chord  give  qualitatively  accurate 
results.  The  spanwise  pressure  distribution  atX  =  2/3  of  the  root  chord  is  given  in  Figure  44 
and  reflects  the  same  trends. 

Figures  45-49  refer  to  a  mode  of  a  half  sine  wave  in  the  chordwise  direction.  The  conver¬ 
gence  of  the  generalized  force  is  shown  in  Figure  45  and  is  seen  to  be  far  better  than  the 
moment  due  to  rotation  (Figure  39)  or  the  lift  due  to  plunge  (Figure  34).  Figure  46  presents 
the  velocity  potential  distribution  near  the  root  chord  and  shows  excellent  convergence  for 
large  numbers  of  boxes  and  fairly  good  results  for  as  few  as  twenty  boxes.  Figure  47  compares 
the  velocity  potential  results  for  twenty  and  fifty  boxes  at  midspan  and  again  shows  good  agre- 
ment.  The  pressure  distributions  near  the  root  and  at  midspan  for  twenty  and  fifty  boxes 
are  shown  in  Figure  48  and  it  also  shows  that  using  as  few  as  twenty  boxes  probably  gives 
reasonable  agreement  with  the  results  for  greater  numbers  of  boxes.  Figure  49  gives  the 
spanwise  pressure  distributions  at  two  different  x  locations  and  shows  that  the  results  for 
twenty  and  fifty  boxes  are  in  fair  agreement,  but  that  using  as  few  as  twenty  boxes  causes 
so  me  discrepancies  near  the  leading  edges. 

Figures  50-53  refer  to  a  mode  shape  of  a  full  sine  wave  in  the  chordwise  direction.  The 
generalized  force  results  in  Figure  50  appear  to  give  excellent  convergence.  The  velocity 
potential  and  pressure  results  in  Figures  51,  52,  and  53  give  results  similar  to  those  for  the 
one  half  sine  wave,  except  that  the  spanwise  pressure  plots  in  Figure  53  show  better  agreement 
at  x  =  50  percent  of  the  root  chord  rather  than  75  percent,  as  was  the  case  for  the  half  sine 
wave. 
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SECTION  V 

COMPARISON  WITH  EXPERIMENT 


As  noted  earlier  in  this  report  the  input  and  output  data  of  the  Sonic  Box  program  can  be 
interpreted  in  either  of  two  ways,  both  ultimately  giving  the  correct  results. 

.  f  Apj  Zi  dS 

(1)  Pressure  coefficients  p/q)  and  generalized  force  coefficients  (  J  — - )  can 

be  taken  directly  from  the  output  data  if  the  input  data  is  understood  to  be  s 


2 

(a)  -j—  as  a  function  of  %  $  for  pressures 

«  A  as  a  function  of  x,  y  for  generalized  forces. 

(2)  If  the  input  data  is  understood  to  be  z  as  a  function  of  x,  y  then  both  the  output  for 
“pressures”  and  “generalized  forces”  must  be  divided  by  the  root  chord  b  to  give  proper 
magnitudes  for  Ap/q  and  f  Apj  f  \  dS  respectively. 

Ps  q~S 

The  second  interpretation  in  this  section  will  be  used  denoting  the  input  data  z  as  a  function 
of  5f  and  y.  to  compare  predicted  results  with  experimental  data  for  the  aspect  ratio  3  rec¬ 
tangle  of  Figure  11,  the  aspect  ratio 2delta  of  Figure  27,  and  the  cropped  deltas  of  Figures  30 
and  31. 


1.  THE  ASPECT  RATIO  3  RECTANGLE 

Figure  11  presents  the  model  used  for  unsteady  pressure  measurements  reported  in 
Reference  3.  Figure  12  presents  the  first  bending  mode  shape  which  was  excited  at  26  CPS 
and  used  for  the  measurement  of  unsteady  pressures  at  the  root,  50,  70,  and  90  percent  of  the 
span.  The  pressure  measurements  were  actually  done  over  a  Mach  number  range  from  0.24  to 
1.30,  but  can  be  compared  with  the  Sonic  Box  program  only  for  Mach  1.0. 


Reference  3  presents 


a  spanwise  plot  of  the  section  lift  coefficient  Cg  ,  defined  by: 

ah 


( Pg  ~  Pu  )  dy 
qb  ah 


(  Pj>  -  P  u  )  d  x 
q  b  (  uf  h  f  /  V  ) 


q  K  hf 

In  the  bending  mode,  the  generalized  force  coefficient,  L,n,  is  defined  to  be: 

f  /b(R£  -  pu  > 1  <y> d  *  d?  /  *  (7)  [/  ( h  ~pu]  dl0  d* 

0  0  x  u _  _  _0 _ _ _ 

L"  =  ps  qS 


(6) 


(7) 


In  terms  of  C*  this  is: 
£ah 


L 


M 


(8) 
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Now,  for  a  given  mode  shape  z(y),  the  Sonic  Box  program  prints  out  a  generalized  force 
coefficient  which  is  equal  to  b  L^.  Therefore: 

‘-:  =  b'-,,=1rLA<»>c.  <V=^T-f  *,>c.  <15-  (91 


bkhT  r~ 

~s  /  *' 


/  z(y>c»  dy  = 


dy 


Using  the  mode  shape  of  Figure  12  for  z*(y)  and  the  data  of  Reference  3  for  C«  (Jr)  a  numerical 

s  ah 

evaluation  of  the  integral  gives  /  ?(y)Cfl  dy  =  4.0  ft.  Then  the  experimental  value  of 

*0  £flh 


L*  i8  4-Ok  h  =  (4.0)(.244)^ 
^11  18  s  1  (2.29) 


=  1.75.  Figure  54  plots  this  number  against  the  various 


values  of  L*  ^  printed  out  by  the  program  for  varying  number  of  boxes  along  the  root  chord. 
The  predicted  value  is  roughly  50  percent  higher  than  the  experimental  value. 

Reference  3  gives  experimental  pressures  in  the  form: 


qcu  hy/V 


(Pl  ~  Pu  >  b 
qk  hT 


Therefore: 


n  Ap*Ap 

Again,  for  a  given  mode  shape  z(y)  the  Sonic  Box  program  prints  out  ( -  J— )  =  b  { — —■)- khTCL 

M  M  1  *7] 


Since  k  =  0.244  and  h.  =  4.098  was  selected,  (Ap/q)*  =  C  was  obtained.  Therefore  the 

1  P«h 

quantities  printed  out  as  “pressure”  are  directly  related  (by  a  convenient  choice  of  h^}  to 
Cn  .  Figure  55  is  a  plot  of  the  theoretical  and  experimental  chordwise  pressure  distributions 


of  midspan.  The  predicted  values  are  of  the  correct  order  of  magnitude,  but  are  considerably 
in  error  over  the  forward  part  of  the  chord.  Figure  56  shows  a  similar  comparison  for  other 
spanwise  stations. 

The  discrepancies  between  theory  and  experiment  appear  to  be  the  result  of  two  major 
factors  (aside  from  thickness  and  boundary  larger  effects) : 

a.  The  linear  differential  equation  is  not  strictly  applicable  for  rectangular  wings,  but 
instead  should  be  applied  to  low  aspect  ratio,  highly  swept  wings,  at  “high  reduced  frequen¬ 
cies.” 

b.  The  basic  assumption  in  the  Sonic  Box  method  of  constant  <f>  within  any  one  box  is  not 

applicable  at  the  leading  edge  where  —  ®.  Figures  13-15  illustrate  how  this  error  causes 

the  predicted  values  of  £  to  be  exaggerated  near  the  leading  edge.  When  these  values  are  in 
turn  “smoothed”  by  a  least  squares  process,  the  initial  error  causes  a  waviness  which  can 
distort  the  slope  of  <£  over  the  chord.  This  effect  is  particularly  noticeable  near  10-30  percent 
chord  and  near  the  trailing  edge. 

2.  ASPECT  RATIO  2.0  DELTA 

Figure  27  is  a  sketch  of  an  aspect  ratio  2.0  delta  wing  which  was  the  subject  of  unsteady 
pressure  measurements  reported  in  Reference  4.  The  model  was  oscillated  in  rotation  about 
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an  axis  at  approximately  71  percent  of  the  root  chord.  Pressure  measurements  were  taken 
along  a  line  at  36  percent  of  the  semispan  on  one  surface  only. 

Reference  3  presents  the  data  in  the  form: 

|Cp|  =  I  : p™°".  I  cm: 


In  terms  of  (p^ 


Pu>/9 


it  is  shown  to  be 

1  c  I  -  I  1  _^“pu  l 

1 

1/  Ap  \*l 

1  CP  1  "  2  1  q  I 

2b 

1  \  4  /  1 

(12) 


where  (Ap/q)*  =  b(Ap/q)  is  the  value  of  “pressure”  printed  out  bv  the  Sonic  Box  program. 
The  experimental  results  are  shown  in  Figure  57a  in  the  form  I C  Pc2 1  =  icp/£l ,  a  in  degrees. 


These  values  should  be  represented  by  the  Sonic  Box  program  by 


r  c  i  -  1 

l  Ap  \*  | 

!  pa  1  2  b  a 

\  Q  /  ! 

(13) 


In  the  calculations  for  this  case,  the  modal  polynomial  was  given  by: 

z  =  -  0.006179  +■  0.017446  ”x 


(14) 

(dimensions,  ft) 


Therefore  a  =  =  -0.017446  radians  =  -1°.  Also,  since  in  this  case  b  =  0.5  ft  it  is  shown 

that  2ijjg  =  -1.0  and  the  experimental  values  of  I  Cp^l  will  be  directly  comparable  to  (Ap/p)* 

printed  out  as  “pressure”  by  the  program.  The  comparison  of  predicted  and  experimental 
values  is  shown  in  Figure  57a  and  is  seen  to  be  poor.  Figure  29  shows  a  similar  comparison 
between  the  Sonic  Box  method  and  Jones’  low  aspect  ratio  theory.  The  comparison  between 
the  two  theoretical  methods  is  seen  to  be  quite  good. 


Figure  57b  sheds  some  light  on  the  disagreement  between  theory  and  the  experimental  data 
of  Reference  3.  The  experimental  data  shown  in  Figure  57a  were  taken  by  a  slender  probe 
which  traversed  the  right  hand,  upper  surface  of  the  experimental  model.  Those  results  are 
also  shown  (for  50  percent  chord)  as  the  empty  circles  in  Figure  57b.  One  is  inclined  to  suspect 
considerable  interference  between  the  probe  and  the  model,  particularly  at  Mach  1.0.  Other 
data  were  taken  with  the  probe  on  the  upper  surface;  however,  in  this  case  a  slender  “fuselage” 
section  was  added  to  the  lower  surface  of  the  model.  Those  results  are  shown  as  the  solid 
circles  in  Figure  57b.  Finally,  a  third  set  of  data  were  taken  with  orifices  within  the  model 
itself.  Those  results  are  shown  as  the  triangles.  A  prediction  of  the  Sonic  Box  is  also  shown 
in  Figure  57a  for  comparison. 

The  considerable  uncertainty  within  the  experimental  data  leads  one  to  minimize  the  im¬ 
portance  of  the  difference  between  theory  and  experiment  by  questioning  the  accuracy  of  the 
measurements  themselves  while  continuing  to  accept  the  limitations  of  the  analytical  method. 


3.  THE  CROPPED  DELTAS 

Reference  5  reported  the  measurement  of  oscillatory  moments  on  a  family  of  cropped  delta 
wings.  Two  of  the  test  cases,  the  aspect  ratio  1.5  and  2.0  wings,  will  be  used  here  for 
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comparison  with  the  Sonic  Box  program.  The  two  wings  under  consideration  are  shown  in  Fig¬ 
ures  30  and  31.  Reference  5  presents  the  experimental  data  in  the  form: 


Reo I ( M ) 
2qS  ba 

imaq  (M  ) 
2qS  b  a  kb 


_  Imaq  (M)  b 
2  q  S  a  k  (b)2 


bu>  T 


where:  =  — —  ,  b  =  mean  chord,  ui  =  frequency  in  radians/sec.  Now,  for  a  given  rotational 

mode  shape,  z'fx)  =  +  A(0x  ,  the  Sonic  Box  program  can  print  out  a  generalized  force 

coefficient  defined  to  be 

L*=bL»=  Ap(T,  y)  z  ir,7)dS  =  ^  /  Ap  (K.y  KA^+A^Til 


A,0  b 


/  Ap  (*x  ,  y  )  (*x- x o  ) 
JS 


where  is  the  location  of  the  rotational  axis,  given  by  =  -A^/A^.  Finally  if,  in  accor¬ 
dance  with  Reference  5,  the  nose  up  pitching  moment  is  called  M,  it  is  related  to  by: 


-A,0  bM 


Noting  that  a  =  -A1QI  mQ  and  m  ^  are  obtained: 

4-  Real  (L*) 

=  - - - - - 

a 

m.  =  -±-l-mqJ  j L  *  J 
a  2  k  ( b  -i  A^ 

For  instance  for  the  aspect  ratio  1.5  cropped  delta  calculations  input  data  are: 
'  b  =  .5833  ft 
b  =  .3333  ft 
A1q  =  1.832 
k  =  0.0953 


Therefore 


+  Real  (  L„  )  /  1.305 


m.  =  +-  Imaq  CL,,  )  /  0.  07106 
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For  the  aspect  ratio  2.0  cropped  delta  similar  calculations  yield: 

ma  =  Real  (L*)  /  1.305 
*  imag  (L*)  /  0.0649 

where 

b  =  .5104  ft 
b  =  .2917  ft 


2.094 


k  =  0.087 


The  results  are  presented  in  Figures  38-61. 


(23) 

(24) 


Figure  58  is  a  plot  of  the  moment  coefficients  mQ  and  versus  Mach  number.  The  open 

circles  represent  experimental  data  from  Reference  5.  Also  shown  are  the  Sonic  Box  predic¬ 
tions  for  Mach  1.0.  In  Figure  58  the  Sonic  Box  predicts  values  for  ma  and  which  are 

approximately  20  and  70  percent  higher  than  experiment,  respectively,  for  the  aspect  ratio 
1.5  cropped  delta  rotating  about  its  forward  axis.  The  experimental  data  shown  was  generally 
taken  with  an  amplitude  of  2°.  however,  the  dashed  line  for  mQ  represents  data  taken  for 

1°  amplitude.  The  large  discrepancy  between  predicted  and  experimental  data  for  is  typical 

of  other  linearized  analysis  (Reference  5). 

Figure  59  is  a  similar  plot  for  the  aspect  ratio  1.5  cropped  delta  rotating  about  its  aft  axis. 
The  predicted  mQ  coefficient  is  still  roughly  20  percent  higher  whereas,  agreement  on 

has  improved  from  the  forward  axis  case. 


Figure  60  refers  to  the  aspect  ratio  2.0  delta  rotating  about  its  forward  axis.  Again,  the 
mQ  coefficient  is  the  right  order  of  magnitude,  whereas,  is  off  by  100  percent.  The  same 

comments  apply  for  the  aft  rotation  axis,  Figure  61. 
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SECTION  VI 

CONCLUSIONS,  RECOMMENDATIONS 


1.  CONCLUSIONS 

a.  Interpretation  of  Mode  Shape  Information 

The  existing  Sonic  Box  computer  program  contains  normalizing  and  non-dimensionallzing 
factors  which  may  not  be  apparent  to  the  user.  When  mode  shapes  are  fed  into  the  programs 
as  coefficients  A  or  as  deflections  at  discrete  coordinates,  there  are  two  interpretations 

ran 

which  yield  correct  results,  adopting  a  dimensional  coordinate  system  x,  y,  zf. 


(1)  To  obtain  pressure  coefficients  Ap/q  directly  from  the  program,  put  in  modes  in  the 
■jf  r  Ap:  7 j  dS 

form  t-  as  a  function  of  xandy.  To  obtain  generalized  force  coefficients  J  —  qs —  directly 

S  ^ 

from  the  program,  put  in  modal  data  in  the  form  *2/  y/~b  as  a  function  of  xf  and  y. 


(2)  If  modal  data  is  assumed  to  be  in  the  form  lasa  function  of  “2  and  y,  then  both  the 
pressures  and  generalized  forces  must  be  divided  by  the  root  chord  b. 


b.  Smoothing  of  Input  and  Output 

The  Sonic  Box  program  uses  smoothing  procedures  in  two  ways;  first,  to  fit  polynomial 
surface  to  modal  deflection  data  which  has  been  put  in  at  discrete  coordinates  on  the  surface; 
and  second,  to  fit  polynomial  surfaces  to  the  calculated  values  of  velocity  potential. 

In  putting  in  deflection  data,  one  coordinate  point  is  sufficient  for  plunge  motions;  two  points 
are  sufficient  for  rotation  or  roll  motions;  and  five  points  are  sufficient  for  a  chordwise  half 
sine  wave.  In  putting  in  a  full  sine  wave  and  3/2  sine  waves  it  was  found  that  proportionately 
more  coordinate  points  were  required.  That  is,  using  nine  points  on  a  full  sine  wave  and 
thirteen  points  on  a  3/2  sine  wave  did  not  give  the  same  input  accuracy  as  five  points  on  a 
half  sine  wave. 

The  smoothing  technique  for  calculated  velocity  potentials  wo  lies  quite  well.  It  can,  however, 
be  distorted  by  large  local  errors  in  <f>  near  the  leading  edge  which  are  reflected  over  the 
whole  surface. 


c.  Convergence 

For  the  aspect  ratio  1.5  delta  wing  at  least  thirty  boxes  along  the  root  chord  are  required 
to  obtain  convergence  of  lift  due  to  plunge  and  other  generalized  forces.  Even  for  large 
numbers  of  boxes,  there  was  some  oscillation  of  the  pressures  near  the  trailing  edge. 

For  the  aspect  ratio  3  rectangular  wing  in  its  first  bending  mode,  as  few  as  fifteen  boxes 
along  the  chord  gave  converged  pressures,  but  that  the  converged  solutions  were  inaccurate 
at  the  trailing  edge  and  from  10  to  30  percent  of  the  chord. 

For  the  aspect  ratio  2.0  delta  wing  more  than  forty  boxes  along  the  root  chord  would  be 
required  to  obtain  a  converged  solution  for  moment  due  to  rotatioh.  Predicted  pressures  \yere 
in  fair  agreement  with  Jones’  low  aspect  ratio  theory. 
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For  the  aspect  ratio  l.S  and  2.0  cropped  deltas  more  than  forty  boxes  along  the  chord  would 
be  required  for  converged  solutions  for  moment  due  to  rotation.  Since  the  moments  for  aft 
rotation  axes  were  more  easily  converged  than  those  for  forward  axes  the  greatest  uncertainty 
lies  over  the  aft  part  of  the  surface. 

For  the  70-degree  delta  about  thirty  boxes  are  required  to  converge  pressures  and  lift  due 
to  plunge.  For  higher  numbers  of  boxes  the  results  could  not  be  distinguished.  For  moment 
due  to  rotation  about  the  apex,  fifty  boxes  were  not  sufficient  for  convergence.  This  implies 
the  oscillation  of  pressures  near  the  trailing  edge.  In  rotation,  the  velocity  potential  seemed 
to  be  reasonably  well  converged  for  as  few  as  twenty  boxes,  leading  one  to  conclude  that  the 
uncertainty  is  introduced  after  the  calculation  of  potentials.  For  a  mode  shape  of  a  half  sine 
wave  in  the  chordwise  direction  potentials,  pressures  and  generalized  forces  were  reasonably 
well  converged  for  as  few  as  twenty  boxes.  Similar  conclusions  hold  for  the  mode  shape  of  a 
full  sine  wave  in  the  chordwise  direction. 


d.  Comparison  with  Experiment 


For  the  aspect  ratio  3  rectangle  in  its  first  bending  mode,  the  predicted  values  of  the 


generalized  force,  ff  Ap(  x,y  ) 7G, y  )  d^dy* 

S 


was  50  percent  higher  than  the  experimental 


value.  The  predicted  pressures  were  of  the  proper  order  of  magnitude  over  most  of  the  aft 
part  of  the  wing,  but  were  inaccurate  at  the  leading  edge,  trailing  edge,  and  the  forward 
quarter  of  the  surface.  Strictly  speaking,  the  theory  does  not  apply  to  this  unswept  planform. 
While  part  of  the  error  is  inherent  in  the  linearized  analysis,  the  planform  treated,  and 
disregard  of  thickness  and  boundary  layer  effects,  part  can  also  be  attributed  to  inadequacies 
in  the  assumed  form  of  the  velocity  potential  near  the  leading  edge. 


For  the  aspect  ratio  2.0  delta,  agreement  with  experimental  pressure  measurements  was 
poor,  whereas,  agreement  with  Jones’  low  aspect  ratio  theory  was  fairly  good.  There  is  some 
evidence  to  indicate  a  certain  amount  of  uncertainty  within  the  experimental  data. 


For  the  aspect  ratio  1.5  and  2.0  cropped  deltas,  the  in  phase  moment  coefficient,  ma ,  was 

g  enerally  fairly  well  predicted.  The  out  of  phase  coefficient,  ,  was  anywhere  from  20  to 

100  percent  high  in  magnitude.  These  results  were  typical  of  other  linearized  analyses,  for 
instance,  those  in  Reference  5. 


2.  RECOMMENDATIONS 


a.  The  existing  Sonic  Box  program  should  be  modified  to  print  out  pressures  and  general¬ 
ized  force  coefficients  which  eliminate  dimensional  confusion  and  do  not  require  further 
adjustments  by  the  user.  Simple  revisions  to  automatically  divide  the  pressure  and  generalized 
forces  by  the  root  chord  would  require  the  user  to  always  put  in  modal  data  in  the  form  zf  as 
a  function  of  Tt  and  %  where  and  2  are  dimensional  quantities. 

b.  More  complicated  revisions  should  be  undertaken  to  improve  convergence.  These 
would  include  improved  representation  of  the  geometry  of  leading  edge  of  the  wing  (now 
jagged)  and  improvements  in  the  allowable  form  of  the  velocity  potential  in  any  one  box  (now 
assumed  constant). 

c.  Since  there  is  considerable  doubt  about  the  value  of  existing  experimental  data  and 
improved  measurement  and  data  analysis  techniques  are  rapidly  coming  available,  unsteady 
pressure  measurements  should  be  made  on  oscillating  wings  in  transonic  f?ow.  This  would 
result  in  a  fair  test  of  existing  methods  and  provide  impetus  for  improvements. 
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Figure  5.  The  Aspect  Ratio  1.  5  Delta  Wing 
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Figure  7.  Real  Part  of  the  Velocity  Potential  for  an 
Aspect  Ratio  1.  5  Delta  Wing  in  Rotation, 
y  =  1  ft  =  26. 67  Percent  Semispan 
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Figure  8.  Imaginary  Part  of  the  Velocity  Potential  for  an 
Aspect  Ratio  1.  5  Delta  Wing  in  Plunge,  y=l  ft 
=  26.67  Percent  Semispan 
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Figure  14.  Chord  wise  Distribution  of  the  Imaginary  Part 

of  the  Velocity  Potential  in  First  Wing  Bending, 
19,  24,  and  28  Boxes 
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Figure  15.  Chordwise  Distribution  of  the  Imaginary  Part  of  the  Velocity  Potential  in  First 
Wing  Bending,  32  Boxes 
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Figure  16.  Raw  and  Smoothed  Values  of  the  Imaginary  Part  of  the 
Velocity  Potential  in  the  First  Row  of  Boxes 
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Figure  17. 
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Figure  22.  Chordwise  Pressure  Distribution  in  Plunge; 
24„  25,  and  26  Boxes 
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Figure  23.  Chordwise  Pressure  Distribution  in  Plunge;  34,  35,  and  36  Boxes 
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Figure  24.  Chordwise  Pressure  Distribution  in  Plunge;  44,  45,  and  46  Boxes 
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Figure  25.  Chordwise  Distribution  of  the  Imaginary  Part  of  the  Pressure  in  First  Wing 
Bending:  6,  11,  15,  and  19  Boxes 
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Figure  28.  Convergence  of  the  Generalized  Force  in  Rotation 
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Figure  33.  Convergence  of  the  Moment  due  to  Rotation  with 

Increasing  Number  of  Boxes  Along  the  Root  Chord 
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Figure  35.  Smoothed  Values  of  the  Imaginary  Part  of  the  Velocity 
Potential  in  the  First  Row  of  Boxes 
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Figure  36.  Chordwise  Distribution  of  the  Imaginary  Part  of  the 
Pressure  in  the  First  Row  of  Boxes 
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Figure  38.  Spanwlse  Distribution  of  the  Imaginary  Part  of 
the  Pressure,  5?  =  0.667  bD 
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Figure  40.  Smoothed  Values  of  the  Velocity  Potential  in 
the  First  Row  of  Boxes 
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Figure  43.  Chordwise  Distribution  of  the  Real  Part  of  the  Pressure  at  Midspan 
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Figure  44.  Spanwise  Distribution  of  the  Real  Part  of  the 
Pressure  along  5c  =  0.  667  of  the  Root  Chord 
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Figure  47.  Chordwise  Distribution  of  the  Real  Part  of  the  Velocity 
Potential  Along  the  Midspan  for  a  Mode  Shape  of  a  Half 
Sine  Wave 
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Figure  48.  Chordwise  Distribution  of  the  Real  Part  of  the  Pressure 
For  a  Mode  Shape  of  a  Half  Sine  Wave 
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Figure  49.  Spanwise  Distribution  of  the  Real  Part  of  the  Pressure 
for  a  Mode  Shape  of  a  Half  Sine  Wave  in  X 
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Figure  50.  Convergence  of  the  Generalized  Force  with  Increasing 
Number  of  Boxes  Along  the  Root  Chord 
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igure  51.  Chordwise  Distribution  of  the  Real  Part  of  the 

Velocity  Potential  for  a  Mode  Shape  of  a  Chordwise 
Sine  Wave 
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Figure  52.  Chordwise  Distribution  of  the  Real  Part  of  the  Pressure  for 
a  Mode  Shape  of  a  Chordwise  Sine  Wave 
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Figure  53.  Span  wise  Distribution  of  the  Real  Part  of  the  Pressure 
For  a  Mode  Shape  of  a  Full  Sine  Wave  in  x 
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Figure  54.  Absolute  Value  of  the  Generalized  Force  in 
First  Wing  Bending 
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Figure  55.  Comparison  of  Chordwise  Pressure  Distribution 
With  Experiment 
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Figure  57.  Data  for  the  Aspect  Ratio  2  Delta 

(a)  Comparison  of  Theory  and  Experiment 

(b)  Uncertainty  in  Experiment 
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Figure  59.  Comparison  of  the  Predicted  Moment  Coefficients  with  Experiment 
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Figure  61.  Comparison  of  Predicted  Moment  Coefficients  with  Experiment 
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TABLE  m 


INPUT  DATA,  70*  DELTA  WING 


Mode 

Aoo 

A10 

u(cps) 

b(in) 

(in/sec) 

Plunge 

0.02 

0.0 

12.0 

50.0 

13675.0 

Rotation 

0.02 

-0. 0006 

f! 

It 

II 

x(in) 

Defl 

1/2  Sine 

0.0 

0.0 

tf 

ri 

If 

12.  b 

0.01414 

25.0 

0.02 

37.5 

0.01414 

50.0 

0.0 

Sine 

0.0 

0.0 

If 

It 

ft 

6.25 

0.01414 

12.5 

0.02 

18.75 

0.  01414 

25.0 

0.0 

31.25 

-0. 01414 

37.5 

-0.02 

43. 25 

-0. 01414 

50.0 

0.0 

3/2  Sine 

0.0 

0.0 

ft 

IT 

11 

4.167 

0. 01414 

8.333 

0. 02 

12.50 

0. 01414 

16.67 

0.0 

20.83 

-0. 01414 

| 

25.00 

-0. 02 

29.17 

-0. 01414 

33. 33 

0.0 

37.50 

0. 01414 

41.67 

0. 002 

45.83 

0. 01414 

50.0 

0.0 

i 

i 
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TABLE  IV 

INPUT  DATA,  ASPECT  RATIO  1.5  DELTA 


Mode 

Ann 

00 

Aio 

u)  (cps) 

b(ft) 

aoo 

Plunge 

1.0 

0.0 

3.183 

10.0 

1000. 0 

Rotation 

0.0 

0.1 

If 

»f 

Tt 

TABLE  V 

INPUT  DATA,  ASPECT  RATIO  3  RECTANGULAR  WING 


Mode 

mm 

DefI 

Km 

— 

Bending 

0.0 

0.0 

26.5 

mm 

1023.  5 

0.229 

0.143 

0.457 

0.389 

0.616 

0.738 

0.915 

1.15 

1.14 

1.6 

1 

1.37 

2.09 

1 

1.60 

2.58 

II 

1.83 

3.11 

2.06 

3.61 

2.  2866 

4. 098 

■ ■ 

TABLE  VI 

INPUT  DATA,  ASPECT  RATIO  2  DELTA  WING 


r~  —  ~ 

Mode 

A00 

A10 

o><cps) 

b(ft) 

aoo  <ft/sec> 

Rotation 

-0. 00617 

■ 

33.7 

0.5 

i  nnn  n 

*  VWi  V 
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TABLE  Vn 

INPUT  DATA,  CROPPED  DELTAS 


AR 

Mode 

Axis 

> 

o 

o 

A10 

b(ft) 

(ft/sec) 

1.5 

Rotation 

Fwd 

1.832 

0.5833 

1000.0 

1.5 

?! 

Aft 

1.832 

0.5833 

tt 

2.0 

Tf 

Fwd 

2.094 

0.5104 

II 

2.0 

?! 

Aft 

-0.  8125 

2.094 

0.5104 

If 
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APPENDIX  I 

A  METHOD  FOR  INTERPRETATION  OF  MODE  SHAPE  INFORMATION 


Adopt  the  convention 

z  =  deflection  with  units  of  length 
x  =  coordinate  ”  ”  ”  ” 

y  _  H  It  It  II  II 

The  coefficients  A  define  a  polynomial  in  X  and  y  which  in  turn  defines  a  mode  shape. 

But  suppose  it  is  not  known  if  our  deflections  and/or  coordinates  have  been  normalized  by 
dividing  by  the  root  chord  or  some  power  of  the  root  chord,  b.  Then  the  correct  form  of  the 
polynomial  must  be  found. 

The  way  to  proceed  is  to  make  two  successive  runs  with  the  computer  program.  In  each 
run  a  plunge  mode  and  a  rotation  are  put  in  about  the  leading  edge.  Keep  the  modal  coefficients 
constant  between  the  two  cases;  keep  geometric  similarity;  but  change  the  root  chord.  By 
comparing  the  pressures  and  generalized  forces  obtained  in  the  two  runs  the  correct  form  of 
the  polynomial  is  determined. 

1.  GENERALIZED  FORCES 

a.  Method  for  Polynomial  Coefficients 


Suppose  the  coefficients  A^  define  a  polynomial: 

The  powers  r  and  t  are  unknown  but  must  be  found.  Define  mode  1  as  a  plunge  mode  by: 

=A°° 

Define  mode  2  as  a  rotation  about  the  leading  edge  by: 


«o  bt 


Now  evaluate  the  generalized  forces  L21  and  L22,  the  lift  and  moment  due  to  the  rotation 
respectively.  41  44 


Lzi  =  qS  /  P27i  dS  =  qs dS  =  ~^S~  /  p2  dS  =  qs  L 


=  <T_«z 

l22  =  ■kfp*r*dS=isfiAJ~'  7dSs%rf  p2rdS  = 


AIOb 

qS 


=  A^b 


»b  <v=A«> 


r-t+! 


CMaa* 


Preceding  Page  Blank 
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a»  =  .  -s.  =  Ai0  b 


for  zero  frequency 


L2i  =  Aoo  b  CLaAl0b  =AooAio  cLa  b 


r-tfl 


.2  Jdr- 2t4l 


L22=A|0b  CM  A,0  b  -  A|C  b‘ 


Thus,  the  generalized  forces  have  been  found  in  terms  of  the  coefficients  and  A1Q,  the 


root  chord  b.  C  and  CM  and  the  unknown  powers  r  and  t. 

Lj  q  (J 


Two  successive  runs  are  made  on  the  computer  with  zero  frequency.  In  each  run  the  geomet¬ 
ric  similarity  (CT  and  C,.  constant),  is  kept  while  the  coefficients  A^  and  Ain  are  held 

Mq  oo  XU 

constant;  only  the  root  chord  b  is  varied. 

Compare  the  generalized  forces  for  runs  I  and  II. 

L«‘/L®  =  (b,I,/bIIr,)2r-' 


L«>  /L««>  = 


2  r  -  2t  +  I 


Abbreviate  the  above  equations  by: 


r  and  t  are: 


=  2  (' 


2  r  -  f 

RL=  rb 

rm=  Rir"2t  41 

log  Rl  log  Rm 


log  R0  log  R0 


—  ) 
i 


t  =  1  + 


ioq  Rl  log 

log  Rb  log  Rb 


For  example.  Rg  =  10,  R^  =  10,  and  RM  =  100,  then  r  =  0.5  and  t  =  0.  So  the  correct  form 
of  the  modal  polynomial  is : 

2  _  v  .  ~m~n 

J/fe  "  X  Amn  x  y 

b  m,  n 

b.  Method  for  Discrete  Deflections 


Numbers  to  represent  coordinates  and  deflections  are  used,  for  instance 

Deflection:  z/br 
Coordinate:  x/bt 
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The  powers  r  and  t  are  unknown,  that  is,  it  is  unknown  if  the  deflection,  coordinates,  or  both 
are  interpreted  by  the  program  to  be  divided  by  the  root  chord  or  some  power  of  the  root 
chord.  To  answer  these  questions  zero  (or  very  small)  frequency  plunge  and  rotation  modes 
were  put  in  by  using  discrete  deflections  at  two  coordinates  for  each  mode. 

Plunge '.  7, /b  =  S,  of  7  /  b*  =0.0 

=  at  -0 

Rotation'.  72/br  =0.0  at  *x  /  b^  *0.0 

=  -t3 

Now  evaluate  the  generalized  force  coefficients  L91  and  L00,  where  1  refers  to  plunge  and  2 
refers  to  rotation.  ** 

La  =  3s-/  *.  dS  8,  •>'  ds  =  Qj  P*  IS  =^2  .  s,  cLz 


8>  b'  cLa  “2  •  8,  c>-a  ^  8,  “r  cLa  (S2bV^9b,)  =  Ml  cLn  b 


2  r  -t 
t 
a 

r-t  +  l 


* dS = iff  )  *  % - c“ 


s 

*5  .r-t+l 

.  s2b  -  dr 


=  CM„  <V#  f  b 


2r  -2t-H 


0 


’a  dY 

Two  successive  runs  were  made  on  the  computer  keeping  geometric  similarity  (CT  and  C.. 

^Qt  ^  CL 

remain  constant) ;  keeping  the  same  numbers  8  ,  8.  and  0  ;  but  changing  the  root  chord  from 
I  II 

b  to  b  .  The  relationship  between  the  generalized  force  for  the  two  runs  are 


/L 


21 


U)_(n),2r-t 


=  (b  /b  ) 


l«/lS*  blIVb0II)2r'2,+ 


These  are  the  same  equations  used  for  the  polynomial  coefficient  case  and  the  same  equations 
can  be  used  to  solve  r  and  t. 

2.  PRESSURES 

a.  Method  For  Polynomial  Coefficients 
Suppose  the  coefficients  Amn  define  a  polynomial 


T  =  E  A  (xVbM  (y/bM 


m.n 


Define  mode  1  as  a  plunge  mode 


=  A 


oo 
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Define  mode  2  as  a  rotation  about  the  leading  edge 


br  =  A'°  hf 


Now  at  some  point  on  the  surface  evaluate  the  downwash  in  each  mode. 


W  _  _dX_  ,  i  k  z* 

V  "  dx  b 

(*),'  ik  A."'" 

(v)#s  A10  fcf  f  ( 1  +  ik  ^/b  ) 


u  y 

With  the  assumption  that  -g^-was  known  at  this  point,  evaluate  p. 

dp  dp  /w  \ 

p  =  =  (v ) 


In  each  mode  there  is 


»,  =  -IJ  likA°o  0rH  ! 

p  =  —  {a.  br_  f  )  (  14  ik  x /b  ) 
z  da  w 

Now  if  two  successive  runs  are  made  with  the  computer,  keeping  geometric  similarity, 
keeping  the  same  k,  and  keeping  the  same  coefficients  A  and  A.fl  it  is  found  that  the  ratios 
of  the  pressures  between  “the  runs  are:  u 


Mode  1: 

_«>/  in)  . 
P.  '  p  i  - 

(I)  (I) 

Mode  2: 

s 

ii 

Abbreviating  these  equations, 


r  and  t  are: 


RP,  =  rB 

R_  =  Rr‘ 
P2  B 


ioq  Rp, 
log  Rb 

log  Rp| 
log  Rb 


log  RPfe 

log  Rg 
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b.  Method  For  Discrete  Deflections 
Again  plunge  and  rotation  modes  are  defined 


v>r 

=  8, 

of  "x  / 1 

t  =  0. 0 

and  ft 

vb' 

=  0.0 

at  x/b 

*■  ^  0.0 

CM 

<50 

II 

--ft 

such  that 

P  = 

( 

i  k  8  b 

da 

P2  = 

dp 

da 

\  ft 

4) 

Again  two  successive  runs  were  made,  I  and  n.  The  geometric  similarity  is  kept  and  ft,  8 , 
and  82  are  held  constant;  and  only  b  is  changed. 


These  are  the  same  equations  found  in  part  2.a,  therefore 


r 


t 


I  + 

I  -5- 


lOQ  Rp( 
log  R0 
log  Rpt 
log  Rb 


log  Rp2 
log  RB 
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TABLE  VIE 


INPUT  DATA  SHEET  FOR  THE  INTERPRETATION  OF  MODAL  POLYNOMIALS 


1  12 

24 

36 

48 

60 

72 

80 

23 

0. 0001 

100.0 

1000. 0 

27 

5.0 

2.0 

1.0 

30 

0.0 

100.0 

50.0 

plunge 

46 

1.0 

0.0 

-  87 

1.0 

1.0 

1.0 

1.0 

1.0 

rotation 

46 

0.0 

1.0 

memmm 

23 

1000.0 

■ 

27 

1 

1.0 

1 

30 

Irg 

■■ 

5.0 

fSSn; 

plunge 

46 

rotation 

.  46 

0.0 

1.0 

TABLE  IX 

INPUT  DATA  SHEET  FOR  THE  INTERPRETATION  OF  MODAL  DEFLECTIONS 


1 

12 

24 

36 

48 

60 

72 

80 

23 

0. 0001 

100.0 

1000.0 

27 

5.0 

2.0 

1.0 

30 

0.0 

100.0 

50.0 

87 

1.0 

1.0 

1.0 

1.0 

1.0 

98 

2.0 

2.0 

1.0 

plunge 

101 

0.0 

0.0 

1.0 

1.0 

0.0 

1.0 

98 

2.0 

■KEI 

1.0 

Mm 

rotation 

0.0 

0.0 

105 

1.0 

HkOI 

1.0 

23 

10.0 

1000. 0 

27 

5.0 

2.0 

1.0 

30 

0.0 

10.0 

5.0 

98 

2.0 

2.0 

1.0 

plunge 

101 

0.0 

S3 

1.0 

— 

1.0 

0.0 

1.0 

98 

2.0 

2.0 

1.0 

rotation 

101 

0.0 

0.0 

0.0 

105 

1.0 

0.0 

1.0 
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TABLE  XI 

INPUT  DATA  SHEET  FOR  THE  ASPECT  RATIO  3 
RECTANGLE  IN  FIRST  WING  BENDING 


1  12 

24 

36 

48 

60 

72 

80 

23 

26.5 

1.5 

1023.5 

27 

32.0 

1.0 

1.0 

30 

0.0 

0.0 

2.2866 

87 

0.0 

1.0 

1.0 

0.0 

1.0 

Qfi 

11  A 

A  A  |  V 

1.0 

11  ■  G 

101 

0.75 

0.0 

0.0 

105 

0.75 

0.23 

.14 

\ 

109 

0.75 

0.46 

.39 

I 

113 

0.75 

0.69 

.74 

/ 

117 

0.75 

0.92 

1.15 

f  First  b 

ending  mod 

e 

121 

0.75 

1.14 

1.60 

> 

125 

0.75 

1.37 

2.09 

I 

129 

0.75 

1.60 

2.58 

\ 

133 

0.75 

1.83 

3.11 

j 

137 

0.75 

2.06 

3.61 

f 

-  141 

0.75 

2.2866 

4.098 

' 
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TABLE  XH 

INPUT  DATA  SHEET  FOR  THE  ASPECT  RATIO  2.  0  DELTA  WING  IN  ROTATION 


1  12 

24 

36 

46 

60 

72 

80 

23 

33.7. 

0.5 

1000.0 

27 

20„  0 

l.f 

1.0 

30 

0.0 

O.i, 

0.25 

46 

-0.00617 

-0.017446 

—  87 

0.0 

1.0 

0.0 

0.0 

0.0 

TABLE  Xin 

INPUT  DATA  SHEET  FOR  THE  CROPPED  DELTA  WINGS  IN  ROTATION 


1  12 

24 

36 

48 

60 

72 

80 

23 

26.0 

0.5833 

1000.0 

I 

27 

30.0 

2.0 

1.0 

}  AR  1.. 

5,  forward  axis 

30 

0.0 

0.5 

0.25 

) 

46 

-0.4034 

1.832 

-  87 

0.0 

1.0 

0.0 

1.832 

— 

23 

0.5104 

fmmmm | 

■  ■ 

27 

2.0 

■  I 

II  1 

30 

■ 

0.4376 

■ 

~ _ i£_ 

IHBSKvm 

2.094 

mm 

■ 

—  46 

-0.8125 

2.094 

■■■ 
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